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(Dated: February 4, 2008) 

We report on a study of inclusive B~ and B° meson decays to D°X, D°X, D + X, D~ X, DfX, 
D~X, A+X, A~X, based on a sample of 231 million BB events recorded with the BABAR detector 
at the T(4S) resonance. Events are selected by completely reconstructing one B and searching for a 
reconstructed charm particle in the rest of the event. From the measured branching fractions of these 
decays, we infer the number of charm and anti-charm particles per B decay, separately for charged 
and neutral parents. We derive the total charm yield per B~ decay, n~ — 1.202±0.023±0.040to'o2g , 
and per B° decay, n° = 1.193 ± 0.030 ± 0.0341q'o35 where the first uncertainty is statistical, the 
second is systematic, and the third reflects the charm branching-fraction uncertainties. We also 
present the charm momentum distributions measured in the B rest frame. 
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PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



I. INTRODUCTION 

The dominant process for the decay of a b quark is 
b — » cW*~ [1], resulting in a (flavor) correlated c quark 
and a virtual W. In the decay of the W, the production 
of a ud or a cs pair are both Cabibbo-allowed and should 
be approximately equal, the latter being suppressed by a 
phase-space factor. The first process dominates hadronic 
b decays. The second can be easily distinguished as it 
produces a (flavor) anticorrelated c quark. Experimen- 
tally, we investigate correlated and anticorrelated charm 
production through the measurement of the inclusive B- 
decay rates to a limited number of charm hadron species, 
i.e. £>°, D°, D+, D~, D+, D~, A+, A~, ~ c and char- 
monia, because all other charm particles decay into one 
of the previous hadrons. 

The analysis presented here exploits a substantially 
larger data sample than the original BABAR result [2]. 
It also employs a more sophisticated fitting method to 
extract, in a correlated manner, the number of recon- 
structed B mesons and the charm hadron yields, which 
reduces the experimental systematic uncertainty. Other 
measurements [3-7] of these rates are more statistically 
limited and/or do not distinguish between the different 
parent B states. Besides the theoretical interest [8-11], 
the fact that anticorrelated charm particles are a back- 
ground for many studies also motivates a more precise 
measurement of their production rates in B decays. 

Most of the charged and neutral D mesons produced 
in B decays come from correlated production B — * DX. 
However, a significant number of B — > DX decays are 
expected through 6 — * ccs transitions, such as B — > 
£)(*)|)(*)^(*)( n7r ). Although the_ branching fractions of 
the 3-body decays B — ► D^D^K have been mea- 
sured [12, 13], they do not saturate B — > DX transi- 
tions [2]. It is therefore important to improve the preci- 
sion on the B — > DX branching fraction. 

By contrast, anticorrelated D~ production, 
B — > DjD(nir), is expected to dominate B decays to D s 
mesons, since correlated production needs an extra ss 
pair created from the vacuum to give B — > Df K~ (rnr). 
There is no prior published measurement for correlated 
-D+ production. 

Correlated A+ are produced in decays like 
B — > A+pTr~(ir), while anticorrelated A~ should 
originate predominantly from B — > S C A~ (tt) . The decay 
B — > S C A~ has recently been observed [14], confirming 



*Also at Laboratoire de Physique Corpusculairc, Clermont- 
Ferrand, France 
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the hypothesis of associated S c A~ production. An- 
other possibility for anticorrelated A~ production is 
B — > A+A~K, the baryonic analogue of the DDK decay. 

This analysis uses T(45) — > BB events in which either 
a B + or a B° meson (hereafter denoted £? r0 c'd) decays 
into a hadronic final state and is fully reconstructed. We 
then reconstruct D, D s and A+ from the decay products 
of the recoiling B~ {B°) meson and compare the flavor 
of the charm hadron with that of the reconstructed B 
(taking into account B°-B° mixing). This allows sepa- 
rate measurements of the B~_ (B°) — > D° X. D+ X_, D+ 
X, A+ X and B~ (B°) -> L>° X, D~ X, D~ X, A~ X 
branching fractions. 

We then compute the average number of correlated 
(anticorrelated) charm particles per B~ decay, N~ 



N- = £S(B- 

c 

Nz = J2B(B- 



CX), 

Ox), 



(1) 

(2) 



where the sum is performed over C = 
{D\ D+, D+,__ A+,_ S D , (cc)}_ or C = 
{D°, D , D s , A c , (cc)}, where (cc) refers to all 
charmonium states collectively. We neglect anticorre- 
lated E c production, as it requires both a cs and an 
ss pair in the decay to give Et c f2 c . We then sum 
and iVjT to obtain the average number of charm plus 
anti-charm quarks per B~ decay, n~ = N~ + . We 
similarly define N®, N® and n° c for B° decays. 

The above method also lends itself to a measurement 
of the momentum distribution of each charm species di- 
rectly in the rest frame of the parent meson, because the 
four-momentum of each recoiling B is fully determined 
from those of the T(4S) and of the reconstructed B. The 
resulting charm spectra can then be compared to theoret- 
ical predictions in the same frame [15]. This avoids the 
significant smearing due to the Lorentz boost from the 
parent-_B frame to the T(4S) frame affecting earlier mea- 
surements, such as those reported in [3]. These spectra 
might also show indications of four-quark states [16]. 



II. BABAR DETECTOR AND DATA SAMPLE 

The measurements presented here are based on a sam- 
ple of 231 million BB pairs (210 fb _1 ) recorded at the 
T(45) resonance with the BABAR detector at the PEP- 
II asymmetric-energy B factory at SLAC. The BABAR 
detector is described in detail elsewhere [17]. Charged- 
particle trajectories are measured by a 5-layer double- 
sided silicon vertex tracker and a 40-layer drift cham- 
ber, both operating in a 1.5-T solenoidal magnetic field. 
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Chargcd-particlc identification is provided by the av- 
erage energy loss (dE/dx) in the tracking devices and 
by an internally reflecting ring-imaging Chcrenkov de- 
tector. Photons are detected by a CsI(Tl) electromag- 
netic calorimeter. We use Monte Carlo simulations of 
the BABAR detector based on GEANT4 [18] to optimize 
selection criteria and determine selection efficiencies. 



III. B MESON RECONSTRUCTION 

We reconstruct B + and B° decays (-B re c'd) m the 
modes B+ -> Z)M°7r+, D^°p+, Z)M a+ and 5° -> 
D^-tt+, D < -* ) ~p+, £)W- a +. D° candidates are re- 
constructed in the K + ir 



A' H 



and 
are 
D* 
and 



K°ir + Tr~ (Kg —* tt + it~ ) decay channels, while D~ 
reconstructed in the K + ir~ir~ and K®tt~ modes, 
candidates are reconstructed in the D*~ — > D°tt~ 
D*° -> D°n° decay modes. 

The kinematic selection of fully reconstructed B decays 
relies on two variables. The first is AE = Eg — y/s/2, 
where E B is the energy of the reconstructed B candi- 
date in the e + e~ center-of-mass frame and y/s is the 
invariant mass of the initial e + e~ system. The sec- 
ond is the beam-energy substituted mass, defined by 
TO ES = V( s / 2 + P» -VbYIE'I -p B , where p B is the 
-Brec'd momentum and (£Jj,pj) is the four-momentum of 
the initial e + e~ system, both measured in the laboratory 
frame. We require \AE\ < tictae, using the resolution 
(Jae measured for each decay mode, with n = 2 or 3 
depending on the decay mode. If an event contains sev- 
eral B + (B°) candidates, only the highest-purity B-decay 
mode is retained. The purity is defined, for each i?-decay 
mode separately, as the fraction of signal B decays with 
toes > 5.27 GcV/c 2 , normalized to the total number of 
reconstructed B + (B°) candidates in same interval. 

The signal yield Nb of reconstructed B mesons is ex- 
tracted from a fit to the toes spectra (Fig. 1). The B sig- 
nal is modeled by a Crystal Ball signal function Tcb [19] 
which is a Gaussian peaking at the B meson mass mod- 
ified by an exponential low-mass tail that accounts for 
photon energy loss. The B combinatorial background is 
modeled using the empirical ARGUS phase-space thresh- 
old function T arg [20]. All the signal and background 
parameters in these functions are extracted from the 
data. The signal yields of reconstructed B + and B° 
mesons are N B + = 200359±705 and N B o = 110735±424, 
where the errors reflect the statistical uncertainty in the 
number of combinatorial background events. These num- 
bers provide the normalization for all the branching frac- 
tions reported below. 

The contamination of misreconstructed B° events in 
the B + signal (and vice-versa) induces a background 
which peaks near the B mass. From the Monte Carlo 
simulation, the fraction of B° events in the reconstructed 
B + signal sample is found to be Co = 0.038 ± 0.009(syst), 
and the fraction of B + events in the reconstructed B° 
signal sample c+ = 0.028 ± 0.007(syst). The system- 
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FIG. 1: mES spectra of reconstructed (a) B + and (b) B° 
candidates. The solid curve is the sum of the fitted signal 
and background whereas the dashed curve is the background 
component only. 



atic uncertainties take into account possible differences 
in reconstructing real or simulated events, as well as 
branching-fraction uncertainties for those B decay modes 
contributing to the wrong-charge contamination. 



IV. INCLUSIVE CHARM BRANCHING 
FRACTIONS 



We now turn to the analysis of inclusive D, D, D~, 
Df , i4+ and A~ production in the decays of the B mesons 
that recoil against the reconstructed B. Charm particles 
C are distinguished from anti-charm particles C. They 
are reconstructed from charged tracks that do not belong 
to the reconstructed B. The decay modes considered are 
listed in Table I along with their branching fractions. 
Those are taken from Ref. [21] except in the case of the 
Df — > 4>tt + channel [22] for which we use the more precise 
measurement reported in Ref. [23] . 
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FIG. 2: Charm (left) and anti-charm (right) mass spectra 
in the recoil of B + candidates, for the subsample of events 
with mES > 5.270 GeV/c 2 (B signal region). The solid curve 
shows the result of the two-dimensional fit. The dark shaded 
areas show the contribution of reconstructed D,D, and 
signal in the recoil of combinatorial B^ cc , d background. 
The light shaded area corresponds to the fitted combinatorial 
(anti-) charm background. 



A. Charm particle yields 

The numbers of charm (anti-charm) particles are ex- 
tracted from an unbinned maximum likelihood fit to 
the two-dimensional distribution [m-Eg,m c rg\], where 
mEs is the beam-energy substituted mass of the re- 
constructed B and m c ,£\ is the mass of the charm 
(anti-charm) particle found among the recoil products. 



FIG. 3: Charm (left) and anti-charm (right) mass spectra as 
for Fig. 2 but in the recoil of B° candidates. 



Figs. 2 to 5 show the results of these fits, projected 
onto the m c axis, for events in the mES signal region 
(toes > 5.270 GeV/c 2 ). The probability density function 
used to fit the [m-ES, m c (c)\ distributions is the sum of 
four components : 

• Pgllg ■ reconstructed charm (anti-charm) signal in 
the recoil of reconstructed B signal, 

• Psbkg '■ reconstructed charm (anti-charm) signal in 
the recoil of combinatorial B background, 

• Pggig ■ combinatorial charm (anti-charm) back- 
ground in the recoil of reconstructed B signal, 
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FIG. 4: (left) and DJ~ (right) mass spectra in the re- 

coil of B + candidates, for the subsample of events with 
m,ES > 5.270 GeV/c 2 (B signal region). The solid curve shows 
the result of the two-dimensional fit. The dark shaded areas 
show the contribution of reconstructed Dj , DJ signal in the 
recoil of combinatorial -B+ c , d background. The light shaded 
area corresponds to the fitted combinatorial (anti-) charm 
background. The Gaussian peak at the D + mass accounts 
for reconstructed D + signal [24]. 



• P 



Cbkg 
Bbkg 



combinatorial charm (anti-charm) back- 
ground in the recoil of combinatorial B background, 

These four components are modeled as follows : 

P Bsig( m vs,mc)= r C B{m E s)'xps(rnc) , 
P Bbk 9 g ( m ^s,mc)= T ARG (m ES xps(m c ) , 



- P Bst 9 ( m ES,?'n-c)= TcB(mEs)xpcomb{m>c) , 

P Bbkg( m VS,mc)= r ARGUES )>< Pcomb(m C ) ■ 



(3) 



The function Tcb with all its parameters fixed from 
the fit detailed in Sec. Ill is used to model the recon- 
structed B signal. The combinatorial B background is 
described as in Sec. Ill by an ARGUS function Targ 
whose shape parameter is floated in the fit to allow for 
a possible charm decay-mode dependence of this back- 
ground. A Gaussian function ps(jn c ^) describes the 



FIG. 5: and D s mass spectra as for Fig. 4 but in the 
recoil of B° candidates. 



mass shape of the reconstructed charm signal. Its mean 
is taken from the data. Its resolution, as measured in 
the data, is consistent with that in the simulation and is 
fixed. The combinatorial charm-background distribution 
is fitted with a linear function p C omb( m c (c)) (except for 
the D° — > K ~ 7r + 7r~ 7r + for which a quadratic depen- 
dence is assumed) [24]. 



TABLE I: Charm particle decay modes and branching frac- 
tions. 



B(C - /) (%) 



D 



=; K-n+ 

D° -> K-TT + TV-TT + 
D + -> K-TT + TV+ 

D+ -> 0tt + (0 -> K+K~) 
D+ -> K*°K+{K*° -> K--K+) 
D+ ^K° S K+(K° 3 ^K+Tr-) 
A+ -> pK-% + 



3.80 ± 0.09 
7.48 ±0.31 
9.1 ±0.7 

4.81 ±0.64 (49.3 ±1.0%) 
4.57 ±0.69 (66.51 ±0.01%) 
2.43 ± 0.42 (68.95 ± 0.14%) 
5.0 ± 1.3 
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TABLE II: p*-averaged reconstruction efficiencies ec for each 
charm final state. The errors reflect the limited Monte Carlo 
statistics. 



C - 
~W 

D° 
D+ 
Df 
Dt 
Dt 
At 



f 



K~ 
K~ 
■ K~ 



ec (%) 



» K*° K 
* K° s 



50.2 ± 0.3 
20.1 ±0.2 
33.7 ±0.2 
33.0 ±0.8 

18.0 ±0.5 

31.1 ±0.8 
26.7 ±0.9 



The reconstruction efficiencies for each charm final 
state C — » / (Table II) are computed from the simu- 
lation as a function of p* , the charm-particle momentum 
in the B rest frame, and applied event-by-event to ob- 
tain the efficiency-corrected charm and anti-charm signal 
yields. These arc denoted respectively by N~(C — > /) 
(N°(C -► /)) and N~(C -> /) (N°{C -> /)) and are 
listed in Table III. We then determine the charm and 
anti-charm fractional production rates B c ^ and B c , 



defined as 



B, 



-(o) _ 



N-(°\C^ f)/[N B+(B0) xB{C->f)], 
N-(°\C -+ f)/[N B+ r B0) x B(C -+ /)] , 



(4) 



where N B + (N B o ) is the number of reconstructed B + 
(B°) mesons, and B(C — > f) is the C — > / branching 
fraction reported in Table I. B~ , B~ , £>° and B% arc 
listed in Table III. 



B. Correlated and anticorrelated charm branching 
fractions 



For charged B, the branching fractions for correlated 
and anticorrelated C production are given by : 



B(B- 
B(B~ 



CX) = B- - c Q B\ , 
OX) = B~ - cqB% . 



(5) 



The correlated (anticorrelated) B~ — > CX branching 
fraction is equal to the charm (anti-charm) fractional 
production rate B~ (JB~ ) in the recoil of reconstructed 
B + mesons modified by a small correction term coS? 
(coB®) that accounts for the B° contamination in the 
reconstructed B + sample. The factors B® and B\ depend 
on the measured B° — ► CX and B° — > CX branching 
fractions, and on the B°B° mixing parameter Xd [21]. 
Doubly Cabibbo-suppressed D° decays (D° — > K + ir~ 
and D° — > K + ir + ir~ir~) are also taken into account. 
We combine the results from the different D and D s 
decay modes to extract the final branching fractions 
listed in Table IV. The probability of the correlated 
.Dj" production observed in B~ decays to be due to a 



background fluctuation is less than 5 x 10 4 . 

For neutral B, charm and anti-charm production in the 
recoil of reconstructed B° mesons have to be corrected for 
B°B° mixing to obtain the correlated and anticorrelated 
charm branching fractions : 



B(B°^CX) = ®c ~ Xd(B + B°) _ C+B + ^ 



B(B° -► CX) 



Xd 

B° c - Xd (B° c +B°) 



(G) 



l-2xd 



c+B, 



+ "2 ' 



The correction factors c+B^ and C+B2 account for 
B + contamination in the B° sample and depend on 
the B~ — > CX and B + — > CX branching fractions. 
Combining the different D° and D s modes, we obtain 
the final branching fractions listed in Table IV. 

We also compute the fraction of anticorrelated charm 
production in B decays : 



w(C) 



B(B -► CX) 



B(B -> CX) ± B(B -> CX) 



(7) 



Here, many systematic uncertainties cancel out (track- 
ing, K identification, D branching fractions, B count- 
ing). The results are given in Table V. 

The main systematic uncertainties are associated with 
the track-finding efficiency, the models used to describe 
the tties and m c ^ distributions, and the particle iden- 
tification efficiency For example, the 2.7% absolute 
systematic uncertainty on B(B~ — > D°X) reflects the 
quadratic sum of 1.3% attributed to the track-finding ef- 
ficiency, 1.6% to the description of the m^s distribution 
by the Targ and Tc B functions, 0.8% to the descrip- 
tion of the m c ^ signal distribution by the ps function, 
1.4% to the particle identification, 0.5% to the Monte 
Carlo statistics, 0.3% to Co, and 0.1% to B®. 

The uncertainty affecting the track-finding efficiency is 
estimated with two different methods. The first uses a 
large inclusive sample of tracks with a minimum number 
of hits in the silicon vertex detector. The second relies on 
an e + e~ — > t + t~ control sample. From these, we derive 
a relative systematic uncertainty of 0.8% per track. 

The modeling of the toes distribution by the Tc B and 
the Targ functions affects both the charm signal yields 
and the numbers of reconstructed B mesons used in nor- 
malizing the branching fractions. The corresponding un- 
certainty is dominated by the dependence of the Targ 
shape parameter on the lower edge of the Toes fit range. 
Varying the latter from 5.195 to 5.225 GeV/c 2 yields a 
variation in the branching fraction that is taken as sys- 
tematic uncertainty. This range was chosen such that the 
branching fractions measured in the simulation change by 
±1 standard deviation. 

The uncertainty associated with the description of the 
charm signal mass shape by the ps function translates 
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TABLE III: Charm and anti-charm efficiency-corrected signal yields and fractional production rates. The uncertainties are 
statistical only. 



/ n^f^oir m /"I /"^ 
\_j LLcCcLy IHULLfc! 


C in recoil of B + ,, 

rec a 

N-(C^f) B-(%) 


C in recoil of B + ,, 

rcc a 

N-(C^f) B-(%) 


111 IfcCUll 

N°(C^f) 


nf R° 
B°(%) 


C in recoil of B° ac , d 
N°(C^f) B°(%) 




5898±126 


77.5±1.6 


691±52 


9.1±0.7 


1731±70 


41.1±1.7 


669±44 


15.9±1.0 




11010±383 


73.4±2.6 


1378±214 


9.2±1.4 


3418±239 


41.2±2.9 


1065±159 


12.8±1.9 




1970±131 


10.8±0.7 


513±89 


2.8±0.5 


3044±122 


30.2±1.2 


869±74 


8.6±0.7 


Dt-><t>n + 


85±24 


1.8±0.5 


385±42 


8.1±0.9 


97±21 


3.7±0.8 


227±30 


8.7±1.2 


-^K*°K+ 


78±39 


1.3±0.6 


567±72 


9.3±1.2 


78±28 


2.3±0.8 


306±50 


9.1±1.5 




0±16 


0.0±0.5 


212±39 


6.6±1.2 


48±19 


2.7±1.1 


148±29 


8.3±1.6 


At -^ p R--k + 


288±52 


2.9±0.5 


210±45 


2.1±0.5 


240±41 


4.3±0.7 


124±30 


2.2±0.5 



TABLE IV: B branching fractions. The first uncertainty is statistical, the second is systematic, and the third reflects charm 
branching- fraction uncertainties [21, 23]. 

Correlated Anticorrelated 

C B(B-^CX)(%) B(B°~>CX)(%) B(B-^>CX)(%) B(B° -> CX)(%) 

D° 78.6 ± 1.6±2.7l^ 47.4 ± 2.0 ± 1.5l^ 8.6 ± 0.6 ± CUj^'j 8.1 ± 1.4 ± 0.5t°i 

D + 9.9±0.8±0.5t°:? 36.9 ± 1.6 ± 1-4±|;| 2.5 ± 0.5 ± O.ll^ 2.3 ± 1.1 ± 0.3lg;? 

< 3.9 at 90% CL 

D+ 1.1±H ± 0.1 jig;? 1.5 ± 0.8 ± 0.1±g;| 7.9 ± 0.6 ± 0.4±^ 10.3 ± 1.2 ± 0A±{- 7 3 

< 2.6 at 90% CL 

At 2.8 ± 0.5 ± 0.3tJ;g 5.0 ± 1.0 ± 0.5t^ 2.1 ± 0.5 ± 0.2±g;| 1.6 ± 0.9 ± 0.2±g;| 

< 3.1 at 90% CL 



TABLE V: Fraction of anticorrelated charm as defined in 
Eq. (7). 



Mode 



B decays 



B° decays 



D°X 0.098 ± 0.007 ± 0.001 

D'X 0.204 ± 0.035 ± 0.001 

D~X 0.884 ± 0.038 ± 0.002 

A~X 0.427 ±0.071 ±0.001 



0.146 ± 0.022 ± 0.006 
0.058 ± 0.028 ± 0.006 

< 0.098 at 90% CL 
0.879 ± 0.066 ± 0.005 

> 0.791 at 90% CL 

0.2431^2? ±0.003 

< 0.403' at 90% CL 



difference between these two efficiencies is then taken as 
an estimate of the corresponding the systematic uncer- 
tainty (1.7% relative uncertainty per kaon and 1.3% per 
proton). 

The statistical and systematic uncertainties in Ta- 
ble IV and Table V are computed separately for each 
charm decay mode; correlated errors are taken into ac- 
count when averaging over D° and D s final states. 



C. Average charm production in B decays 



into an uncertainty on the charm reconstruction effi- 
ciency. It is estimated by fitting the simulated charm 
signal with a double instead of a single Gaussian. 

The systematic uncertainties affecting the proton and 
charged kaon particle-identification efficiency are esti- 
mated using D° — > K~ir + and A° — > pir~ samples re- 
coiling against reconstructed B + and B° mesons. The 
D° or yl° signal yields are extracted in a manner similar 
to that described in Sec. IV A, both with and without ap- 
plying the proton or kaon particle-identification require- 
ments. The ratio of these yields on real and simulated 
samples is proportional to the particle-identification effi- 
ciency in the data and the simulation, respectively The 



To extract N c from the results of Table IV, we still 
need to evaluate the B — > S C X and B — > (cc)X branch- 
ing fractions. Because there exists no absolute mea- 
surement of the .E^-decay branching fraction, the ab- 
solute rates for correlated S c production in B decays 
arc unknown [14, 25]. Therefore, following the discus- 
sion in Sec. I, we assume that B(B — > S C X) = B(B — > 
A~X)-B{B -> A+A-K^k)) [26]._Aj-ecent measure- 
ment [27] indicates that B — ► AtA~K decays have a 
branching fraction of the order of 7 x 10~ 4 , and thus can 
be neglected by comparison to Nc~^° (sec also [2]). Wc 
take B(B -> (cc)X) = (2.3 ± 0.3)% [28, 29] and, using 
Eqs. (1) and (2), we obtain for charm production in B~ 
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decays: 



N~ = 0.968 ±0.019 ± 0.0321^, 
iVf = 0.234 ±0.012 ±0.008;°°^, 
n- = 1.202 ± 0.023 ±0.040±g;g||. 



and in B° decays 



A c ° = 0.947 ± 0.030 ± 0.028±g:gl|, 



0.246 ±0.024 ±0.0091^14, 
1.193 ± 0.030 ± 0.034±g;g||. 



The results reported here are consistent [30] with, and 
supersede those of Ref. [2]. The three-fold increase in 
integrated luminosity accounts for the substantial re- 
duction in statistical error. The experimental system- 
atic uncertainties have been similarly reduced, primar- 
ily through the use of the two-dimensional [m-ESi m c (C)) 
fit, which takes correctly into account the correlation be- 
tween the fitted number of reconstructed B mesons and 
the corresponding charm yield. 



D. Isospin analysis 

The main source of anticorrclated D mesons produced 
in B decays is b — * ccs transitions. In these processes 
isospin should be conserved, leading to the expectation 
that : r(B-_-> D°X) = T(B° -> D~_X) and Y{B~ -> 
D~X) = T(B° -> D°X). However, D mesons can also 
arise from D* mesons, whose decay does not conserve 
isospin since the D*° — ► D~tt + channel is kinematically 
forbidden. Thus isospin invariance actually requires : 



Fdir(B 


-^D°X) 


= r dir (B° 


-> D~X) 


Tdir{.B~ 


-> D~X) 


= r dir (B° 


-> D°X) 


T(B~ - 


• D*°X) 


= r(B° -> 


D*~X) 


T(B~ - 


> D*-X) 


= T{B° 


D*°X) 



(8) 



where Tdir{B — > DX) refers to the partial width of B- 
meson decays to D mesons where the D state is not 
reached through a D* cascade decay. Eqs. (8) lead to 
the following relations involving the measured anticorre- 
lated D branching fractions in Table IV : 



B(B~ 
B{B° 



D°X) - B(B a 
Jb 



D~X) 



r B + 



D°X)- 



and : 

x + x* - 

+B{B° 



\[B{B- 



D Q X)^± 



D a X)- 
^B{B° 



B{B- - 

B(B- - 
■* D-X' 



T~B° 

D-X) 



D-X) 



T „0 



(9) 
(10) 

(11) 



where t^/t^ is the ratio of the B + to the B° lifetime, 



x* = B(B- -> D*° + D*-X) [31]. That both Eqs. (9) 
and (10) must be satisfied is a consequence of isospin 
invariance. From these two equations, we extract x* with 
a chi-squared method, and using in addition Eq. (11) we 
calculate : 



B(B~ - 

B d ir{B~ 



D* 



D { 



D*-X) 
\-D~X) 



B dlr {B -> D° + D~X) 
B(B — > D*° + D*~X) 



= 9.1 ±1.5 ±0.6% 
= 2.1 ±1.7 ±0.7% 

< 4.5% at 90% CL 



0.23±g;i| ± 0.09 



< 0.60 at 90% CL 



B(D* 



D^n- 



B d iv{B- -> D° + D-X) and 



Here the first uncertainty is statistical, the second is 
systematic and includes charm branching-fraction uncer- 
tainties, as well as those affecting the values of t^/t^ 
and B(D*~ -> D°tt-) taken from Ref. [21]. The X 2 of 
the fit to Eqs. (9) and (10) is 0.01 for 1 degree of freedom. 



V. CHARM MOMENTUM DISTRIBUTIONS IN 
THE B REST FRAME 

As the four-momentum of the recoiling B is fully deter- 
mined, each reconstructed charm hadron can be boosted 
into the rest frame of its parent B, yielding the p* distri- 
bution of the corresponding (anti-charm) charm species 
in the B frame. The number of C (C) candidates, their 
fractional production rates and the B — > C (C)X branch- 
ing fractions are then determined in each p* bin by the 
same methods as in Sec. IV, separately for B- and B° 
decays. The systematic uncertainties are assumed to be 
independent of p* , except for the error associated with 
the B° (B + ) contamination in the B + (B°) sample : the 
latter is computed bin-by-bin with a relative uncertainty 
on c + and cq increased to 100%. 

Figs. 6 and 7 show the result for correlated and anti- 
correlated D°, D + , D s and A£ production in B- and B° 
decays, respectively. The numerical values are tabulated 
in the Appendix. 

Correlated D° and D + (Figs. 6a, c and 7a, c) are pro- 
duced in several types of transitions : b — > dt~v, b — > cud 
and b — > ccs which explains the fairly large spread of 
their momentum. High-p* correlated -D's are produced 
in two-body decays such as B- — > D°ir- while low mo- 
mentum D's might come from higher multiplicity final 
states such as B — > DDK(Xught) where Xn g ht is any 
number of pions and/or photons. The latter processes 
are also the main source of anticorrclated D° and D~ 
production (Figs. 6b, d and 7b, d) which explains why 
anticorrelated D spectra are softer than their correlated 
counterparts. 

Anticorrelated D~ spectra (Figs. 6f and 7f) have a 
very different shape compared to anticorrelated D spec- 
tra. They are peaked at high p* values which is suggestive 
of the two-body decays B -> and B -> D^D*~. 
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FIG. 6: Momentum spectra, in the B rest frame, of corre- 
lated (left) and anticorrelated (right) charm particles : D° / D° 
(a)(b), D ± (c)(d), Df (e)(f), Af (g)(h). The error bars are 
statistical only. The histogram in frame (f) represents the 
contribution of B~ —* Z)'*' Z)s*' two-body decays assuming 
the branching fractions of Ref. [21] and [23]. 



FIG. 7: Momentum spectra, in the B° rest frame, of corre- 
lated (left) and anticorrelated (right) charm particles : D°/D° 
(a)(b), D ± (c)(d), Df (e)(f), Af (g)(h). The error bars are 
statistical only. 



VI. CONCLUSIONS 



These decays represent a large fraction of the total anti- 
correlated D~ production as shown in Fig. 6. In contrast, 
the corresponding two-body processes B — ► D^D" and 
are Cabibbo-suppressed. 

In the case of anticorrelated A~ production associ- 
ated with S c production, for decays such as B — > 
S C A~ (Xught), the anticorrelated A~ spectra should have 
a cut-off at p* < 1.15 GeV/c. This is actually observed 
in the data, both in B~ (Fig. 6h) and in B° (Fig. 7h) 
decays. 



We have measured the branching fractions for inclu- 
sive decays of B mesons to flavor-tagged D, D s and Af, 
separately for B~ and B°. We observe a significant pro- 
duction of anticorrelated D° and D + mesons in B de- 
cays, with the branching fractions reported in Table IV. 
These results are consistent with and supersede our pre- 
vious measurement [2]. We find evidence for correlated 
Df production in B~ decays, a process which has not 
been previously reported. 

The sum of all correlated charm branching fractions, 
N c , is compatible with 1, for charged as well as for neutral 
B mesons. The numbers of charm particles per B~ decay 
(n- = 1.202 ± 0.023 ± 0-040±g;g||) and per B n decay 
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(n° = 1.193 ± 0.030 ± 0.0341°;°^) arc consistent with 
previous measurements [2, 4, 28] and with theoretical 
expectations [8-11]. 

Assuming isospin conservation in the b — > ccs transi- 
tion, we show that anticorrelated D mesons are mainly 
produced by cascade decays B — ■> D*X — ► DX. 

Finally, the technique developed for this analysis 
allows us to measure the inclusive momentum spectra 
of flavor-tagged D, D s and in the rest frame of the 
B parent, separately in B~ and B° decays, eventually 
providing insight into S-decay mechanisms. 
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APPENDIX : CHARM p* SPECTRA 

This appendix tabulates the measured p* dependence 
of the branching fractions displayed in Figs. 6 and 7. 
In Tables VI to XIII, the first uncertainty is statistical, 
the second is systematic and includes charm branching- 
fraction uncertainties. Within each table, the statistical 
uncertainties are uncorrelated whereas the systematic er- 
rors are fully correlated. 
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TABLE VI: Correlated and anticorrelated D production in B decays. 

correlated prod. anticorrelated prod. 
p* range ( GeV/c) B(B~ -> X C X) (%) g(5~ --> X-X) (%) 
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TABLE VII: Correlated and anticorrelated D + production in B decays. 

correlated prod. anticorrelated prod. 
p* range (GeV/c) B(B~ -> X C X) (%) g(5~ -> X^X) (%) 
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TABLE VIII: Correlated and anticorrelated D s production in B decays. 

correlated prod. anticorrelated prod. 
p* range ( GeV/c) B(B~ -> X C X) (%) -> X-X) (%) 
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TABLE IX: Correlated and anticorrelated At production in B decays. 

correlated prod. anticorrelated prod. 
p* range ( GeV/c) -> X C X) (%) B(B~ -> X^X) (%) 
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TABLE X: Correlated and anticorrelated D° production in B° decays. 

correlated prod. anticorrelated prod. 
p* range ( GeV/c) B(B~ -> X C X) (%) B(B~ -> X^X) (%) 
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TABLE XI: Correlated and anticorrelated D + production in B° decays. 

correlated prod. anticorrelated prod. 
p* range ( GeV/c) B(B~ -> X C X) (%) B(B~ -* X-X) (%) 
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TABLE XII: Correlated and anticorrelated D s production in B° decays. 

correlated prod. anticorrelated prod. 
p* range ( GeV/c) B(B~ -> X C X) (%) B(B~ -> X-X) (%) 
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TABLE XIII: Correlated and anticorrelated A~£ production in B° decays. 

correlated prod. anticorrelated prod. 
p* range ( GeV/c) B{B~ -> X C X) (%) B(B~ -> X^X) (%) 
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